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PREFACE 


This  report  documents  a  study  performed  by  Systems  Research  Laboratories, 
Inc.  (SRL)t  Dayton,  Ohio,  for  the  Air  Force  Aerospace  Medical  Research 
Laboratory  (AFAMRL) ,  Human  Engineering  Division,  Optical  Countermeasure 
program.  This  work  was  performed  under  Contract  F33615-79-C-0500.  The 
Contract  Monitor  was  Mr.  Donald  McKechnie  and  the  Program  Manager  was 
Maj.  Allan  M.  Dickson.  The  SRL  Project  Manager  was  Mr.  Kaile  Bishop. 

The  author  extends  his  appreciation  to  Dr.  Carroll  Day,  Messrs.  Walter 
Summers  and  Maris  Vikmanis  of  the  Human  Engineering  Division  of  the  Air 
Force  Aerospace  Medical  Research  Laboratory,  Wright-Patterson  AFB,  for 
their  valuable  comments  and  discussions. 
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Section  I 
INTRODUCTION 


The  modeling  of  human  performance  in  an  antiaircraft  artillery  (AAA)  system 
has  been  extensively  studied  by  many  investigators  in  the  past  decade  [e.g., 
Kleinman  and  Perkins  (1974),  Phatek  et  al.  (1976),  Kou  et  al.  (1978)].  Most 
of  these  works  dealt  with  the  modeling  of  human  response  in  a  simple  tracking 
task.  In  the  event  of  interrupted  observations,  the  operator's  tracking 
performance  degrades  significantly  during  the  interruption  period  and  poses 
an  appealing  modeling  problem.  The  author  tackled  this  problem  by  degrading 
•several  observer  and  controller  gains  in  the  model  and  proved  to  be  rather 
successful  (Yu  et  al.,  1980).  Efforts  were  then  directed  to  study  human 
response  in  a  manual  tracking  and  firing  task.  In  this  task,  the  operator 
(gunner)  directly  controls  the  gun  turret  and  fires  tracer  rounds  continu¬ 
ously  toward  the  target.  The  gunner  perceives  the  tracer  ending  position 
and  continuously  adjusts  weapon  pointing  in  azimuth  and  elevation  to 
minimize  the  tracer-to-target  error.  In  this  system  mode,  the  gunner  has  to 
play  both  the  role  of  a  tracker  and  a  lead  angle  computer.  The  conventional 
tracking  task  is  greatly  complicated  by  the  inclusion  of  lead  angle  estima¬ 
tion.  Wei  (1981)  developed  an  observer  gunner  model  which  treated  the 
tracer  information  as  delayed  measurements.  The  intent  of  this  paper  is  to 
extend  the  author's  previous  work  to  consider  an  even  more  general  tracking 
and  firing  scenario,  i.e.  to  consider  a  tracking  and  firing  task  subject  to 
external  measurement  interruptions. 

The  interruptions  occur,  in  the  real  world,  through  various  electronic/ 
optical  countermeasures,  weather,  or  terrain  conditions.  In  this  study, 
extensive  manned-simulation  experiments  were  conducted  at  the  Air  Force 
Aerospace  Medical  Research  Laboratory  of  Wright-Patterson  AFB,  Ohio.  A 
typical  helicopter  operational  trajectory  was  used  in  the  experiment.  The 
trajectory  consists  of  three  phases.  During  the  first  phase,  the  target  is 
standing  still  at  certain  altitude  and  is  half  masked  by  some  terrain  con¬ 
figuration.  At  the  onset  of  the  second  phase,  the  target  pops  up  for  a  full 
unmask  flight  and  moves  horizontally.  Blanking  of  target  is  administered  in 
this  phase  only.  Blanking  durations  range  from  1.5  sec,  3.0  sec,  6.0  sec, 


and  full  blanking.  Certain  repetition  of  blanking  duration  is  also 
Included. 

The  structure  of  the  gunner  model  in  a  tracer-directed  fire  system  in  Wei 
(1981)  is  adopted  and  generalized  here.  Nonlinear  ballistic  equation  is 
used  to  compute  the  loci  of  elevation  projectiles.  The  model  consists  of  a 
reduced-order  observer,  a  linear  feedback  controller,  and  a  noise  remnant 
element.  The  remnant  function  which  lumps  all  of  the  random  effects  due  to 
measurement  noise  and  human  neuromotor  response  noise  is  assumed  to  be 
Gaussian  with  its  covariance  being  a  function  of  estimated  target  velocity 
and  acceleration. 

The  effect  of  observation  interruption  is  modeled  by  exponentially  degrading 
the  observer  gain,  the  controller  gains  pertained  to  observed  states,  and 
the  bias  term  in  covariance  function.  Model  parameters  and  time  constants 
are  identified  separately  with  respect  to  no-blanking  and  blanking  empirical 
data  via  a  least-squares  minimization  algorithm.  The  computer  simulation  of 
the  designed  model  shows  that  the  model  predicted  tracking  and  tracer  errors 
are  in  good  agreement  with  empirical  data  over  various  blanking  conditions. 


Section  II 

AAA  TRACER-DIRECTED  FIRE  SYSTEM 


In  a  tracer-directed  fire  mode,  the  gunner  perceives  both  the  tracking 
error  and  the  tracer  error  on  a  two-dimensional  visual  display.  The  tracking 
error  is  the  difference  between  the  target  angle  6^  and  the  barrel 
pointing  angle  9g.  It  is  also  referred  to  as  "lag  angle"  later  in  this 
report.  The  tracer  error  e^  is  the  difference  between  the  target  angle  0^, 
and  the  projectile  ending  angle  8  .  In  the  simulation  experiment,  the  pro¬ 
jectile  flight  path  ended  at  the  range  of  the  target.  Each  tracer  round 
disappeared  at  this  point,  0^,  from  the  display.  The  detailed  description 
of  the  configuration  is  described  in  Wei  (1981).  We  will  briefly  summarize 
the  underlying  dynamic  system  in  this  report.  Figure  1  is  the  block  diagram 
of  an  AAA  tracer-directed  fire  system. 
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Figure  1.  Block  Diagram  of  an  AAA  Tracer-Directed  Fire  System 


At  any  given  time,  the  target  trajectory  input  eT  is  fed  into  a  visual 
display  device  and  combined  with  the  barrel  pointing  angle  0g,  as  well  as 
the  projectile  ending  angle  0  ,  to  form  error  signals  e^  and  e^.  The  human 
operator  observes  these  error  signals  and  generates  a  control  output  u  via  a 
controller,  or  H-grip,  displacement.  The  control  signal  then  drives  the 
barrel  and  rate  control  plant  for  a  new  barrel  pointing  angle  0  .  Tracer 
round  is  fired  at  this  angle  and  passes  through  the  projectile  ballistics 
computation  to  obtain  the  projectile  ending  angle  0  .  The  task  of  the 
gunner  is  to  constantly  align  the  projectile  ending  angle  to  the  target 
angle,  i.e.  to  minimize  the  tracer  error  The  dynamics  for  the  elevation 

and  the  azimuth  firing  system  are  very  similar.  In  addition,  the  elevation 
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system  can  be  decoupled  from  the  azimuth  system.  However,  the  azimuth 
system  cannot  be  separated  from  the  elevation  system  due  to  a  coupling 
factor  cos  (Qg)EE  the  measurement  equation. 

By  introducing  a  state  vector  x^t)  =  [x^Ct),  xi2(t),  xi3<t)]r,  "T"  means 

"the  transpose  of,"  with  x^Ct)  =  9 iT < c )  “  xi2(t)  *  0iT(t)  "  6iP(t) 

and  x^3  ■  B^Ct),  i  =  1»  2*,  the  following  system  and  measurement  equations 

which  represent  the  underlying  tracer-directed  fire  system  can  be  derived, 
see  Wei  (1981) . 

+  u1(t)  +  E^t)  ui(t-t)  +  Ft  0iT(t)  +  G^t)  (1) 

and 

^(t)  =  C^t)  x^t)  i  =  1,2  (2) 

where 


*If  not  otherwise  specified,  the  first  subscript  index  i  represents  the 
elevation  (i  *  1)  or  azimuth  axis  (i  =  2),  while  the  second  index  repre¬ 
sents  the  i-th  element  or  row  of  a  matrix. 


c2  *  cos  01B(t) 

e1(t)  *  -1.34  x  (1-t)  x  [1+  (0.0052t+0.000486t2)  sin  ©1B(t-T)J 
e2(t)  =  -1.28  x  (1-t) 

g^(t)  ■  (0. 0052+0. 000972t)  x  t  x  cos  0^g(t-T) 
g2(t)  =  0 


;T’  ui’  ^ii*  and  yi2  <*enote  c^e  elevation  or  azimuth  components  of  the 
target  acceleration,  the  gunner's  control  output  and  the  observed  tracking 
error  (lag  angle)  and  tracer  error,  respectively. 


If  we  introduce  a  transformation  on  the  states  x^  and  x^2  by  x^  =  c^x^, 
x^2  *  c^  xi2>  xj^  =  then  Equations  (l)-(2)  can  be  rewritten  as  follows. 


A|(t)  x^  +  u±(t)  +  E^(t)  Uj  (t-r) 


+  F, 


01T(t) 


+  G' (t) 


^(t)  -  C'  x’(t) 


( 


where 


and  x^(t)  »  jxj^t),  x^(t)»  *^3 (t)J Notice  that  the  coupling  factor 
cos  01B  is  removed  from  the  measurement  equation  and  absorbed  into  the 
system  equation.  Equations  (3) -(4)  represent  a  nonhomogeneous  linear 
time-varying  system  with  a  time-varying  delay  in  the  control. 


Section  III 

AN  AAA  GUNNER  BLANKING  MODEL 


In  Wei  (1981),  the  author  proposed  an  observer  gunner  model  for  gunner 
performance  in  a  tracer-directed  fire  system.  The  function  of  the  gunner 
can  be  decomposed  into  two  parts  to  be  modeled.  In  the  first  part,  the 
gunner  observes  continuous  signals  and  makes  an  estimate  of  system  states 
based  on  his  internal  model  of  target  motion.  In  the  second  part,  the 
gunner  utilizes  the  observed  and  estimated  states  to  form  and  exercise  a 
control  action  in  order  to  achieve  his  objective.  The  former  one  corre¬ 
sponds  to  an  estimation  process,  while  the  latter  corresponds  to  a  control 
process.  The  reduced-order  observer  is  used  in  conjunction  with  a  linear 
feedback  control  law  to  model  the  gunner's  function.  The  structure  of  the 
model  is  shown  in  Figure  2. 


Figure  2.  Block  Diagram  of  an  AAA  Gunner  Model 

This  model  structure  is  retained  for  the  blanking  case  except  the  ballistic 
equation  is  no  longer  parameterized  by  a  linear  equation  relating  0p  and  0fi. 
Instead,  a  more  realistic  nonlinear  ballistic  equation  is  used  as  shown  in 
Equations  (3)-(4).  In  addition,  time-varying  gains  are  used  to  model  the 
effect  of  observation  interruption.  We  will  discuss  the  no-blanking  case 
first,  then  the  blanking  case. 


NO  BLANKING 


The  equation  representing  the  gunner's  internal  model  of  the  tracking  and 
firing  system  can  be  written  as 

i^(t)  -  A^(t)x^(t)  +  B^(t)ui(t)  +  E^(t)  u1(t-x)  +  Gj(t)  (5) 

Zi(t)  -  x|(t)  i  -  1,2  (6) 

Since  both  xj^  and  are  measurable,  the  only  state  that  needs  to  be 
estimated  is  x^.  The  state  reconstructor  equation  for  can  be  derived 
by  applying  the  reduced-order  observer  theory  (Luenberger,  1971) 


*  •  • 

Xi3(t)  -  - ^11+^i2) cixi3 (t>  +  kilyil(t)  +  ki2yi2(t) 

-bikilCiUi(t)  "  biki2ciei(t)u1(t~T)  -  ki2ci8i(t) 

*  kilCiCi"lyil(t)  -  kilCiCi_lyi2(t)  (7) 

The  objective  of  the  gunner  Is  to  minimize  the  tracer  error  so  that  a  maximum 
probability  of  hit  could  result.  In  other  words,  the  gunner's  response  in 
the  control  process  would  be  to  stabilize  the  underlying  system,  especially 
the  tracer  error  x^(t);  therefore,  a  linear  feedback  control  law  of  the 
following  form  is  designed  to  achieve  this  objective. 

UfU)  m  Li  (t)  +  vA(t)  (8) 

where 


[Y 


il’  Ti2’ 


Yi3 1 


is  a  vector  of  controller  gains  to  be  identified, 


x^(t)  -  [y^Ct),  yi2(t),  xi3(t)]T 


Is  a  vector  of  measurable  states  and  estimated  state,  v^(t)  is  a  remnant 
noise  function  assumed  to  be  Gaussian  with  zero  mean  and  a  covariance 
function 


r*  a  A  " 

[vi  (t)  V1  (8)]  *  11  +  a12  |eiT  (t)|  +  a13  |9lT(t)|  6(t-8) 


for  all  t  and  s.  are  nonnegative  model  parameters  to  be  determined. 

0iT  and  0iT  are  estimated  target  angle  rate  and  acceleration,  respectively. 
Equations  (7)  and  (8)  represent  the  gunner’s  response  in  the  estimation  and 
control  process  of  the  tracking  and  firing  task.  If  we  define  a  new  state 
vector 

X^,(t)  ■  ^12^  *  xi3^  *  xi3^  ~  *13^^ 

then  the  state  equation  of  the  closed-loop  system  is  obtained  by  combining 
Equations  (7)  and  (8)  with  Equations  (3)  and  (4)  of  the  actual  tracking  and 
firing  system. 

X^t)  -  A^OX^t)  +  D^t)  5^  (t-r)  +  (t)vi(t) 

+  E11(t)vi(t-x)  +  R^t)  (K 


where 


v1  w  «  .  ***■  :V.  ;> 


There  are  seven  model  parameters  in  total,  i.e.,  k^  Yil»  Yi2»  1^3  ai2» 

and  that  need  to  be  determined  from  the  empirical  no-blanking  tracking 
data. 

BLANKING 

The  optical  display  of  target  was  blanked  periodically  according  to  the 
duty  cycles  and  durations  listed  in  Table  1. 


TABLE  1.  BLANKING  CONDITIONS 


Condition 

Duty  Cycle 

(X) 

Blanking  duration 
(sec) 

1 

25 

1.5 

2 

25 

3.0 

3 

25 

6.0 

4 

50 

1.5 

5 

50 

3.0 

6 

50 

6.0 

7 

75  * 

1.5 

8 

75 

3.0 

9 

75 

6.0 

10 

100 

1.5 
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The  duty  cycle  is  defined  as  the  ratio  of  the  blanking  duration  to  the  cycle 
time.  The  blanking  duration  is  the  length  of  time  that  the  target  is 
blanked  so  that  the  subject  cannot  see  the  target.  The  blanking  always 
occurs  at  the  last  portion  of  a  cycle  and  may  reoccur  periodically  over  the 
entire  TOW  firing  period.  An  example  of  blanking  sequence  is  given  in 
Figure  3. 

BLANKING 

PERIODS 


4.0  8.0  12.0  16.0  20.0  T(SFC) 

Figure  3.  Sequence  of  Blanking  for  Condition  4 

The  gunner's  performance  deteriorates  considerably  under  observation 
interruption  via  blanking  the  target.  In  Yu  (1981),  the  effect  of  blanking 
on  the  gunner's  tracking  performance  was  modeled  successfully  by  degrading 
the  gunner's  estimation  gain  k(t)  and  the  controller  gain  y(t).  A  similar 
approach  is  adopted  here  to  model  the  effect  of  blanking  in  a  more  complex 
firing  task.  More  specifically,  the  observer  gain  k^  and  controller  gains 
Yjj  and  y^  which  pertain  to  the  observed  states  xj^  and  x^  are  assumed  to 
decrease  exponentially  as  the  blanking  starts  and  to  increase  exponentially 
as  the  blanking  stops  (see  Figure  4) . 
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During  the  recovery  period:* 


k±(t)  “  ki^cl^  + 

,ki(t0)  ‘  kl(tl)] 

/  t-t.  \ 

ail(t)  “  ail(tl> 

1  -  exp  -  J 

\  /  J 

1  - 


(15) 


(16) 


The  time  constants  x  associated  with  each  gain  parameter,  an<* 
«^(t^)  are  determined  from  the  empirical  tracking  data  collected  in  the 
blanking  experiments,  as  shown  in  the  next  section. 


*In  the  simulation  program,  the  length  of  the  recovery  period  is  defined  as 
the  minimum  of  1.5  sec  and  one-third  of  blanking  period  to  avoid  covariance 
being  negative. 


Section  IV 

PARAMETER  IDENTIFICATION  AND  SIMULATION 


The  least-squares  identification  program  developed  in  Wei  (1981)  was 
modified  to  identify  the  no-blanking  parameters.  Equation  (10)  can  be  first 
decoupled  and  then  approximated,  via  the  Average  Approximation  Method,  by 
the  following  ordinary  differential  equation  (Banks  and  Burns,  1978). 


W±(t)  -  NjU)  W^t)  +  M±  ^(t) 


xi3(t)  “  eix^t^ 


x14(t)  -  -k±(t)c1  x14(t)  +  e1T(t) 


(17) 

(18) 
(19) 


where 


cici  +  biciYil(t) 


biCiYi2(t) 


-1 

CiCi 


CiYil(t)ei  CiYi2(t)ei 


1 

0 

1 

T 

T 

0 

1 

0 

T 

1  0 

’xil(t) 

0  1 

Xi2(t) 

0  0 

Si  * 

x^1(t-x) 

0  0  . 

X'^U-T) 

1 

T 
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(i  +  biYi3)  c±  xi3(t)  -  biCiY13xi4(t)  +  b1c1vi(t) 


% 


CiXi3(t) 


Ciei(t) 


Y13xi3(t-T)  - 


Y13xi4(t-T) 


Vi(t-T)| 


Ci8i(t> 


IDENTIFICATION  OF  MODEL  PARAMETERS 


The  equacion  which  governs  the  mean  of  states  is  obtained  by  taking 
expectation  of  Equation  (17): 


W^t)  -  N±(t)  W^t)  +  M£  i±(t) 


(20) 


where 


wA(t) 


(1+biYi3)cixi3(t)  -  b1ciYi3xi4 


= 


c1*i3(t) 


ciei(t)Yi3 


Xi3(t'T) 


Xi4(t-T)} 


ci8i(t) 


The  first  and  second  component  of  represent  the  model  prediction  of 
ensembled  mean  of  tracking  and  tracer  error,  respectively.  On  the  other 
hand,  the  covariance  matrix  P^(t)  satisfies  the  following  equation: 


P^t)  -  N^t)  P^t)  +  P^t)  N^U)  +  L^t)  £±(t)  L^t) 

where 

^(t)  -  E 


(21) 


21 


! 


bici 


0 

0 


ciei(t) 


0 

0 


Vt} 


ail(t)  +  ai2 


- 

- 

V‘> 

+  ai3 

6iT(t) 

- 

0 

ail(t"T)  +  ai2 

9iT(t-T) 

+  ai3 

9  iT(t_T) 


The  first  and  second  diagonal  element  and  pi22(t)  of  Zj  (O  represent 

the  square  of  the  model  prediction  of  standard  deviation  of  tracking  and 
tracer  error,  respectively.  Notice  that  time-varying  parameters  are  assumed 
for  kp  y^t  Yi2,  and  to  reflect  the  effect  of  blanking.  Since  the 
blanking  effect  to  the  estimation  of  target  velocity  xi3(t)  is  predominantly 
expressed  through  degradation  of  k^,  there  is  no  need  to  consider  a  time- 
varying  Yi3f  «12»  and  a^. 


The  steady-state  value  of  the  parameters  are  first  identified  via  a  least- 
squares  curve-fitting  identification  program.  The  reference  curves  to  be 
fitted  are  obtained  from  empirical  tracking  and  tracer  data  collected  in  the 
manned  simulation  experiments  without  observation  interruption.  These 
experiments  were  conducted  on  an  AAA  simulator  at  the  Air  Force  Aerospace 
Medical  Research  Laboratory.  Three  simulated  helicopter  trajectories 
ranged  1500  M,  2000  M,  and  2500  M  from  the  AAA  system  were  used  as  target 
trajectories.  Figure  5  shows  some  characteristics  for  the  1500  M  trajec¬ 
tory.  Let  x^(t)  and  x^2(t)  t*ie  empirical  ensemble  means  of  tracking 
and  tracer  errors  and  s^(t)  and  si2(t)  be  the  corre9Pon<*ing  standard 
deviations.  These  empirical  means  and  standard  deviations  were  obtained  by 
averaging  and  computing  the  variance  of  the  empirical  data  from  40  simula¬ 
tion  runs  with  the  same  target  trajectory  and  the  same  subject. 
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The  parameters  were  identified  by  minimizing  the  cost  function 

defined  as  follows: 


i  =  1,2 


where  tg  is  the  initial  time  when  a  selected  tracer  round  reaches  the  range 
of  the  target,  t^  is  the  time  when  the  last  tracer  round  is  fired,  is  a 
positive  weighting  factor  chosen  to  be  one  in  the  identification  runs. 

The  direct  search  algorithm  developed  in  Wei  (1981)  was  modified  to  identify 
the  steady-state  values  of  the  parameters.  The  tracking  and  tracer  data  of 
the  helicopter  trajectory,  ranged  1500  M  from  the  AAA  simulator,  without 
blanking,  were  used  to  obtain  the  following  steady-state  parameter  values 
shown  in  Table  2. 

The  time  constants  associated  with  the  parameters  were  determined  empirically 
from  the  data  of  the  1500  M  trajectory  with  blanking  condition  5  and  listed 
in  Table  3. 
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TABLE  2.  STEADY-STATE  PARAMETER  VALUES 


*Insensitive  for  elevation  case,  no  degradation  is  necessary. 


Notice  that  the  time  constants  for  t.  ,  t  ,  and  t  are  the  same  for  both 

k  Y2  °1 

elevation  and  azimuth  gunner  model.  This  is  as  expected  because  the  gunner 
manipulates  the  H-grip  for  elevation  and  azimuth  tracking  indiscrlminantly 
with  respect  to  the  observation  interruption.  On  the  other  hand,  a^(tg) 
and  a^(t^)  for  the  azimuth  case  are  considerably  greater  than  that  for  the 
elevation  case.  This  reflects  the  steeper  increase  of  uncertainty  to  the 
target's  position  along  the  azimuth  axis,  because  the  azimuth  component  of 
target  acceleration  is  much  higher  than  the  elevation  component. 

SIMULATION  RESULTS 


The  gunner  model  was  implemented  on  a  CDC  CYBER  175  computer  to  simulate  the 
man- in- the- loop  AAA  tracking  and  firing  task.  For  the  convenience  of  numeri¬ 
cal  computation.  Equations  (20)  and  (21)  are  discretized  into  the  following 
form: 


<j>n  +  h"  £n 
-i  -i  M 


(23) 


,n  n 

4-  P, 


(24) 


where 


r  • 

I- 


H 


■f 


®xp (t)  •  o]  do  •  M. 

l  n  — i 


*:  i 


exp[Ni(tn)  •  a]  4a  •  (tn> 


Sd  -kK  > 


2"  ‘  W 


A  =0.06  seconds 


A  simulation  program  was  developed  which  uses  the  recursive  Equations  (23) 
and  (24)  to  simulate  a  closed-loop  AAA  tracking  and  firing  task.  Inputs  to 
the  simulation  program  are  the  time  history  of  range  and  acceleration  of  the 
target  aircraft,  the  initial  angular  position  and  velocity  of  the  target, 
the  number  of  blanking  intervals,  and  the  blanking  intervals  in  chronological 
order.  Outputs  of  the  simulation  program  are  model  predicted  mean  tracking 
error  and  its  standard  deviation. 

Simulation  results  are  shown  in  Figures  6  through  17  for  the  blanking 
conditions  3,  4,  5,  6,  7,  and  10.  The  solid  curves  in  these  figures  are  the 
empirical  data  which  are  obtained  by  averaging  the  results  of  40  experi¬ 
mental  runs.  The  dashed  curve  is  the  model  prediction  of  ensembled  mean 
and  standard  deviation. 


Figure  6  and  Figure  7  show  the  comparison  of  model  versus  empirical 
elevation  mean  and  standard  deviation,  azimuth  mean  and  standard  deviation 
of  both  tracking  errors  (lag),  and  tracer  errors  for  the  no-blanking 


case.  Figure  10  and  Figure  11  show  the  results  for  blanking  condition  5 
which  had  a  SO  percent,  3.0  second  blanking  occur  during  [11.01,  14.01] 
seconds . 

Of  particular  Interest  is  the  comparison  of  the  empirical  standard  deviation 
in  Figure  6  and  Figure  10.  The  effect  of  blanking  the  target  to  gunner's 
performance  is  clearly  demonstrated  by  the  sharp  increase  of  the  standard 
deviation  of  tracking  errors  during  the  blanking  period  [11.01,  14.01] 
seconds.  This  effect  is  very  well  modeled  by  degrading  selected  model 
parameters  as  Indicated  in  the  model  prediction  curve  in  Figure  8.  Similar 
agreements  between  the  empirical  data  and  the  model  prediction  can  be  found 
in  Figures  8,  9,  and  12  through  17. 

These  figures  show  that  the  designed  gunner  model  can  provide  consistent 
prediction  of  the  gunner's  empirical  tracking  data  as  well  as  the  tracer 
error  data  for  both  no-blanking  and  blanking  cases.  These  figures  also 
demonstrate  that,  for  a  given  AAA  weapon  system,  the  same  set  of  parameter 
values  and  associated  time  constants  can  be  used  to  predict  the  human 
tracking  and  tracer  errors  for  all  simulated  blanking  conditions. 

However,  due  to  the  fact  that  the  helicopter  trajectory  has  very  low 
elevation  axis  maneuvering,  these  parameters  may  only  hold  for  similar 
types  of  low  maneuvering  helicopter  trajectories.  Reidentification  of 
these  parameters  may  be  needed  for  other  high  maneuvering  trajectories. 

The  computer  execution  time  of  the  overall  simulation  for  an  18  second 
helicopter  trajectory  takes  about  5.60  cp  seconds  on  a  CDC  CYBER  175 
computer . 


Elevation — No  Blanking 


Figure  7a.  Mean  and  Standard  Deviation  of  Tracking  Error' 
Azimuth — No  Blanking 


Azimuth — No  Blankipg 


Figure  8a.  Mean  and  Standard  Deviation  of  Tracking  Erri 
Elevation — 1.5  Seconds,  50  Percent  Blanking 


Figure  9b.  Mean  and  Standard  Deviation  of  Tracer  Error- 
Azimuth — 1.5  Seconds,  50  Percent  Blanking 


Figure  10b.  Mean  and  Standard  Deviation  of  Tracer  Error- 
Elevation — 3.0  Seconds,  50  Percent  Blanking 


Figure  l la.  Mean  and  Standard  Deviation  of  Tracking  Error 
Azimuth — 3.0  Seconds,  50  Percent  Blanking 
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Azimuth — 3.0  Seconds,  50  Percent  Blanking 
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Figure  12a.  Mean  and  Standard  Deviation  of  Tracking  Error- 
Elevation — 6.0  Seconds,  50  Percent  Blanking 
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SUBJECT  33  SUBJECT  33 

TRAJECTORY:  HELICOP1  TRAJECTORY ‘  HELICOP1 

CASE  6  CASE  6 


Figure  13b.  Mean  and  Standard  Deviation  of  Tracer  Error1 
Azimuth— -6.0  Seconds,  50  Percent  Blanking 


Figure  15b.  Mean  and  Standard  Deviation  of  Tracer  Error 
Azimuth — 1.5  Seconds,  75  Percent  Blanking 


Figure  16a.  Mean  and  Standard  Deviati 
Elevation — 1.5  Seconds,  1 


Elevation — 1.5  Seconds,  100  Percent  Blanking 
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Section  5 
CONCLUSION 

This  report  summarizes  the  modeling  of  a  gunner's  performance  in  a  complex 
AAA  tracking  and  firing  task  under  pseudorandom  observation  interruptions. 

The  highlight  of  the  task  is  that  the  gunner  fires  tracer  rounds  without  the 
aid  of  radar  and  lead  angle  computer.  Furthermore,  the  gunner's  performance 
is  greatly  hindered  by  observation  interruptions  via  blanking  the  target  on 
an  optical  display.  A  blanking  model  is  designed  which  consists  of  a  reduced- 
order  observer,  a  linear  feedback  controller,  and  a  remnant  element. 

The  gunner's  performance  is  parameterized  by  the  controller  and  estimator 
gains,  in  addition  to  the  covariance  coefficient  of  the  remnant.  The  effect 
of  blanking  is  modeled  by  degrading  these  gains  and  coefficients  as  a  func¬ 
tion  of  blanking  duration.  An  exponential  decay  form  is  assumed  for  these 
parameters.  The  associated  time  constants  are  determined  from  empirical 
data  collected  in  the  blanking  experiment.  A  direct  search  method  is  used 
to  identify  model  parameters  systematically  while  minimizing  the  least- 
squares  error  between  the  model  output  and  the  empirical  data. 

Computer  simulation  of  the  proposed  gunner  model  show  that  the  model 
predictions  are  in  good  agreement  with  empirical  data  for  various  blanking 
patterns  using  a  typical  helicopter  trajectory.  These  results  demonstrate 
that  the  model  can  adequately  describe  the  gunner's  tracking  and  firing 
characteristics  in  an  AAA  weapon  system  subject  to  observation  blanking. 

This  gunner  (blanking)  model  has  been  incorporated  into  the  MTQ  series  of 
POOl/OBS  AAA  engagement  models  and  is  designated  as  program  P001/0BS  3/6B. 

This  composite  program  POOl/OBS  3/6B  can  be  used  in  the  evaluation  of  air¬ 
craft  survivability  and  performing  weapons  effectiveness  studies.  Documen¬ 
tation  of  POOl/OBS  3/6B  is  in  preparation  at  the  Air  Force  Aerospace  Medical 
Research  Laboratory  and  will  be  distributed  separately. 


APPENDIX  A 

LISTING  OF  PARAMETER  IDENTIFICATION 
PROGRAM  ELEVATION  CASE 


1 1  g,  .if  u  JU  g  •  krtucnincii 

COMMENT.  **0M16EL  10, 10*1 76  0Z95  ,CY*  IT  IOENT  EL  PARA  USE  Hill  TAAJ» 


rl  N. 

AITACH.TAPEl, 1WL6HELSUBJJ  J,CY*1,MR*1. 


DIMENSION  ALFHA(7),PSI(7|  «  A  (7)  ,  0(7)  <  ALP  HAN(  7) 

0l^Sf<3N_5UM5(*T;TiT)T7(7,7r,r  (7,7>,XI  (7,7) ,  3  (7,7>i  F(  A)  ?<  6,r» 
COMMON /ARRAY/X  (  1000),S(1000),ELOD(lgOO)  ,TX«  10  00)  ,TS(  10  00)  ,  RANI  1 00  D 


CCMMON/S/C 0 ,  CEL  ,NSTP,  NOIN  ,T  0 ,  IPT,  Y10,  Y2  0 


INTEGER  Q, SUE1T.FR 


Lb*U 

OC  1  1*1,7 


READ',  M,N017,T0,  IPT,  Y100.Y10 

~pRrNr"T^7»ci;  koi  m,  to  ,iPT,na - 

FORMAT  (1H1, "NO. OF  PTS  ■  “  ,  14,  ZX  , “ORDER*  ",12. 

TKT'rarrres  -,Gl?.5/71X7'*RE:«rEVERir%r?7"_'‘0IW~,-,Y10=I  -,512 

C.S) _ 


DEL*9.03*IPT 


K*K1-KT 
MEAtHl,^: 
C  1*1,  XT) 


r,  rZNNiTTJUMZ,  AZSO,  OU  M3, 


READ  (1,  .31  T  .CUNb.AZ,  AZD,  AZOO.El.ELO,  OUH.AZNN,  0UM2,  AZSO.OJM3 


TXCIDOUMi, 

TS(fG»*3U5S 

COMTIMUt 


55 


$6 


CtSTtNUi 

CO*t.  34 

-|T5T»«-riT=T - 

REA 0»* (ALPHA  (It  *  1*1,7) 


OLOJaAJ 


FCTMAT  (“OALFPAa  G12.5,"I"«G18.5«“I“»5  12. 612.  5, “l~ 


•  3*  »  IDU)3|  1  » 

PRINT  4t.AJ,  IT.SUGIT 

_FC7Hirr^'js'“  -rtn  zrs  ,*r'rr twin  0*5*  n  sr-s  m  rm  inms«”,  i  *r 

OO  100  <*1*7 


Mill  |£Kyl:l  I K3 1 


CO  12  L*l*7 


1F((L.G£.  2)  .AND.  (L.LE.4))  GO  TO  12 

ux^atph™ 

CCNTINUi 

CALL  INTGULFHANrm - 

IF  (AJ.GT.OLCJ)  GO  TO  20 


ULJJ>. J 

PRINT  40*  ALPHAN 


U«I*C«I  ♦£<« 


CC  15  H*l,7 


C  NTlNUc 


ttiiimi 


mi  —  lie  >■ 


GC  TO  25 


IF  ( A  ( <)  .L£.  1.51  A  (<l  *0.0 


CC  33  L*l*7 


k  tiinM 


ITVIiri'Bl'V-l  I 


CM  1.  0 


IF  (CK.  (E.0. 0)60  TO  100 


00  32  N*l,7 

~su(imni*xi  u,m**z“ 

iU'12*SVM2*XI  (2,HI**2 


X1*SQ?T(SUH1) 


trAJ-VLllMM. 


X3*X17X2 


F0TM4T  (*1J«  *,G12.57*  ALPHA*(*,G12.5,*l»,;i2.5**l*,G12.5, 

— 2*i»;6tr.5,*l*T512;y,,«»,;i2.5i*-|*>612.5/ -  -  ~ 

3*  XI ( 11*  »,G1 2.5/*  XI (117X1(2)*  *,G12.5> 


MU*C.3 

F2U,M)»0.0__ 

CCiftNU? 

CC  117  _J*  1  ,J 
oo  Hr  <«  i,  7 
x;<  j,  k  >*o,o 


OC  120  t»l, 

SUi3»SUJj»jA6S < a < in > _ 

If  <SUNT.L£.efS)G0  TO  1000 
00  110  M»l,7 


DC  110  J»l,7 
Zd,  J)*KN)*  V(N.J) 


cciriNUi 

OC  140  J*  1 , 7 


OC  140  L*  1,  7 
CC  140  <*J,7 


XI(J,  L)*XI(JfL»4Z(<,U 
cc^tinui 


SL(4»0 .0 
OC  ISO  J*  It  7 


$im*sui4*xi  <ii  j»»*2 
siH4»sa^r  isuf4i 


LO  15f  J« If  7 
V  (It  J)  »<I  (It  J1/SUH4 


KClNr  f 
0*<00NT-1 

o<r  i?r  <ii7o 

OC  ISO  L*  1 1 7 


F2(2,M)*MQ>»V<K,N)  ♦  F2(Q,M) 
170  CCMUnuJ 
Film.) 


cc  no  r*i(7 


ijgegrinifvgra  #«:«■:  *  i :  t  <re  «  a  **  i •?» « 

l‘|  FJ  ' I 


200  SU15 ( Q)  *S UH5  (Q)  »B (K OUNT, HI  2 

- 31  *3URrjJURT(7in - 

OC  21$  M«l,7 


IF  IOIHT.LE.7)  00  TO  159 

TIT^TTfl - 

SU9IT *  0 


«• 

00  250  <«1 »  7 


.  1600  CALL  EXIT 
.  fcKI 

SU9«0UTINc  1  MGCALPHA  ,EJ) 

. - £ieiM«ff/5RiTS?7SEi«p  {T«o  OTTSEwmroo  ncoarto  ooittx  i  1000  j  ,  ts  uoo  or 

.  1,141U10)0),T  h£T  ( 100  01 


i  •€  c  c  m  i  rj?f. 'A  h  1  n  l  K'iijn  l  i  *  rt  i 1  >  1 1  ri  t* 1 1 


DIMENSION  Ml  4>,P<4,  4)  ,  PI  I  4, 4)  ,P2l 4,  4)  ,  ALPHA  1 7 1  ,  A  (4, 6) 
TTBTIoTTTinTFfTrBl>iT11617T»ll*~t7.>7rAI7rri4,4J,n4)  ,CQC<  4,4)  ,0 141  — 
2.  XT  1 10  91)  ,X4  (10  00)  ,  EOHIIO  30)  fEOOH(lOOO) 

taut lihcs.  <Aa,i),Tri)T7i£Aiifi)VEgii»),i£AjNnitDfEei'iTii))' 
u#T»  XTM./ 


77  ,*).*•  wvn 


57 


c 

c  initialization 


SCAL»CO‘*2/0£L 


1'IMQTU*! 


OC  1  J*l*  NO 


litfUIPI 


CCCiJUO. 


F  (I  I  *  0  • 


IlifVIOniMI 


PI1»1)*0«000  0256279 


■  jfcft-jnmun  iM  i^ij 


X  3< 1) »  0 . 0  00 


ECH(  II  *X3  (II 


1F(ARG.  JT.-200.)  S1«EXP(  ARGI 


C  COMPUTE  AND  STORE  STATES  X3  AND  X4 


9b  FORMAT  (i»G12«  51 


X3(XlT»T3TXIO>EOOfKKI  »DE[T~ 
IF(AL»H1(1I  .EQ.O.I  GO  TO  4 


>  ii«f  l  iSu  firi  fKHll‘lli1UtTP^  iri 


GO  T0  3 


C  COMPUTE  ANO  STORE  ESTIMATED  TARGET  VELOCITY  AN3  ACCELERATION 


ECH(«»*X3(KK)-K4(KKI 


M  1*N0*  •  J 


C  START  INTEGRATION  LOOP 


M  (*l»Y  21 

~ISET*fl - 

CC  110  <<«!»  NSTP 


IFI«.Li.NI  GO  TO  100 


IF(ISET.NE.l)  GO  TO  101 

"SPEA«fi(^Tir-XEMrTRKll**2*(Nr21=TX(in»»r - 

S SIX  SORT  (P(  l'l)l-SEHP(l)  l**ZX  SORT  (P  <2 ,211 -TS<1  ll*»  2 


till  Ut>t.  ♦).001«IS.2*TAU(«)  ♦.'.8b*  1AUI  KK)  **  2)*SINITHtl(  KK«M>  > 


AL9*AIPHA(2) 

j^SWTX 

K  !*<<♦! 

RAaNAA/TlU(KK) 

A2*1.-(TAU(K1)  -TAU  ( KK)  I/O  EL 

"  Cd^*C0*ALE> 

A  (1, 11  ■  •C  OR 

A  ll,»)»-C0*ALPHA(3I 
.  A(E,NC1)*A(1,1|*AL1*A2 

A  (?,  N J I* A  ( 1 ,  21*  ALl*Az 

r% 

00  2  1*1, N02 

« 1*  I  ♦  z 

A ( Jit  I ) *RA 

..  2 

A  (11,  Jt)*-RA 

CCNTINUE 

C 

C 

9* 


CR1*C0*ALPHA  {HI  *AL1*A2 
CRA.tfl^iL^^A  (4) 
SCAL1*S;AL»(AL1*A2I*»2 

TcTT2ti  r.iLPwn?r»«  Axr 


IF((NOO(KK,  100l.EQ.0l  .OR.  (KK.EQ.ll)  PRINT  99,  T  ,M  <11  ♦  Y10,  S1E  AN, 
1  U&rwUTu  liSSd  ' 

FORMAT  (5G12.  5) 


C  COMPUTE  MEA'I  TRACKING  ERROR 


OC  110  I*  1,  NCIM 
.  BcTT?6~ i«1,nCIm 

KIMOIII  *EA  (I,J)*4(J) 

.  TT5  rarran  : 

„  110  CCNTIN'JE 

, - miiT  IV-^CRhIT X3TXK I ♦CRP* X1HXRI - - 

F  (»I*K  JtKKI  ♦CRS*(X<»(KK-MI  -X3(KK-MI  1*0.0  01*(1.-AS»I*(5.2*0.9T2* 

. - C^TAUrKK) )  *  C  O  S  ( THETTKK  -  H1I - — - 

.  00  130  1*1,  NCIM 

. - trur  j*i;z - 

C(I)«C(I)*EA1NT(I,JI»F(J) 

.TTB - COTTOR - 

«  N  (II  *D  IIJ 

. - - - - 

.  130  CONTINUE 

T  ■ 

„  c  COMPUTE  ERROR  DUE  TO  MEAN  TRACKING  ERROR 

.15 - - - 

SMSAN*SNEAN*  (N(1)-XEMP(KKII**2»(H(21-TX  (KP)>**2 

.,15 - 

..  C  COMPUTE  COVARIANCE  MATRIX 

."C - 

CCC(1,1I*  (ALPHA  (5)  ♦ALPHA  C  6)»  ABS  (ECH(KK)  )♦ Al’H A ( 7 1  * AJS (EOOH (KKH  I 

.. - r"»?cTL - - - - 

IFt«.GT.MI  CQC(Z,2I*  ( ALP  HA  (5  I  *  ALPH  A  ( 6)  *ABS  <  E  DH  <  KK-1 1 1 

- 1 — ♦  Al  AHA  <  ?» •  AfS  iEOBUTRK^'O  I  I  *  SCAL 1 - - - 

CALL  MULT  (EA  INT,CQ3  ,  NOIM,  N1,P1,  10) 

.  <,AU  mulT (EA.P,noim,ni,R2,1B) 

.  OC  E20  1*1, NCIM 

. - CC'TirTiTTfTtrH - — 

P(t,JI*Jl(I,  J)*P2(I,JI 


em 

illRDUTt NE  MU.TlEtF.L»LltH,NR> 

- Di^rH;ia»rEnririir,T5H5i  — 

CO  10  I*ltL 

. - im - 

.  00  10  K«1,L 

.  rrwrr; 

OC  5  J*lt Ll»  L 

.  t  eh  f  »  te  hp  *rr  3  t*t  mi 

.5  u-ir»t 

„ - KiFnr-Tnin - 

„  H«K>*T'MP 

„*m - g  KznTznv - 

1F(  MR.  ED.  II  RETURN 


xjc~?o~r»rvr 

00  30  <*I  t L 


TTI"*! r. 

1  I*K 

VC  1 5  7*TTTT7T  “ 

Tt1P*TE  1P»G(  J)*EIIII 

IS 

n»ix*x 

kk»(<*i)*l»i 

L2=L-1 

, - L-c  rtrmrcz - 

.  LJ*I»1 

, - DET-rO  3-*CT,C - 

K1*TI-1I»L»J 

. - Kt*l  J  -  II  *m - 

.  3v  H  (Kll  *H(K3I 

. - m - 

SL1»3UTINE  OSCRT(NOIN«  A«  DEL  iEA«  EA  INTi  NT  I 

- DnrrsnN  Ain.EATtr.  etinriii  ,ccefx301 — 

C  SETS  EA*EXP|A*DELI  .  EAI  NT=  I  NT  EGRAL  EA  0  TO  DEL 

, - NaNitTrrrmx - 

.  NMN3IN»NDIM 

„ - imrmrr=i - 

C  CEF  ( NT)  «1, 

« - OTriT"mTirm - 

II*NT-X 

. - ITT  CGEFTII1  fPEVCOEFTTI+ll  /FLOATTI T - 

c  NT  MUST  6c  AT  LEAST  3 

„ - ~  C«.L"  JTACINDIHtEAINT.TUSOEFttTrCOEri 3IT 

CO  60  L=»3,NT 

, - CALr-T»arrcATEAIIir,T«JIff,HN-,EAl  11 - 

..  IFTl.EO.NTlGC  TO  ro 

.  - 60~ CALL  OTTGTNT]TnvEAXNTTEA«mriCOEF(L11 - 

70  OC  AO  1 1*  1.  N  K.  N0ZN1 

• - rA>ixr»EArm*TTO - 

00  CCNTINUi 

.  - EK1 - 

SUO’DJTINE  0 IAO INOIHf  At  3t  Cl? C2I 

„ - onE'irxoinrnmTi - 

NClNl*  N3IM*1 

- NK»10TN*NTHTI - 

.  NUMDU-l 

- tut - 

IFtCl.El.l.O)  GO  TO  10 

. - or  s- j*r»NNt  rom - 

K«i|*'IKt 


00  6  !■  it  K 

E  TlTTaTTn 

Halt»tiUMl 

fc£TU*N 

TfPS - 


,2.46,2,-0.017964,-9.020511 

1,  .OlHil,  .124433,.  42S18, .  22~44 6F-77.1  P  97  5E-J,.  17302E*~3 


APPENDIX  B 

LISTING  OF  PARAMETER  IDENTIFICATION 
PROGRAM  AZIMUTH  CASE 


.  TToTT.T  JOTI,  3975  TO 0 . — 1760295TWEI.  758-39 60 - 

.  M«P(ON>  . 

.  COMMENT C  W.  6~A  Z I C ,  T  D  *C7  6  fl  2  9 5TCT*  TT“TOEN T  AZ  PUTS  USE- RETI- TRAIT* 
.  FT  N. 


•-  i~i "TWL E"R EUST3 EXT 3 7CT»T f  METs'II - 

AT  TACH.TAPE2, TML6  T  RACERSU3J  S3, 3  T*  1,  10*1.76  0295  ,MR*1. 

FR33*AS  0»T  (TAPE  1, TAP £2, INPUT ,OLTPUT) 

~  Cl  TENS  I TN  Al  FhA  17JTPS I  (71  ,A(7T70{7)  ,ALPWNT71 - 

oHENSIlN  SUf»,(b),E  171, 0(7, 7»  ,Z  (7,7)  ,XI  17 , 7  1 , 3 1 7 , 7 ) ,  F  <  6)  ,  •' 

- ccTmN/APcfl>  /-xriooaTrsriooo)7A70ocioooT,ELT(roonj'i?ANcioo 

1»  TX  (1033)  ,  rs  (1000) 

CC1H3N/S/C0,  lEL  «  NST  P»  NDI  Mi  TOTTPTt  Y1 0 1 Y2  0 , AZ30 - 

uCTIOAl  cm, PAR 

I  NTT  3  F  9  T.SUEIT.FR” - - - 

L  G*  0 

- N  FI  F*7 - 

NFA?t*NJAR-l 

CC  1  1*1*  NP  A  F  - - - 

CC  1  J*l,  NPA  F 

v  (r,  j)*<r  <u  j)*o.o 

IF  II. EX.  II  W(I.J)  =  XI  (I,J)*1.0 

- zcrrjr*T.  o- - 

b(!,JJ*J.: 

T  - CCNTIN'JE  "  - - — - - - - - - 

F  R*  1 

- RE13*,  <1,  NOT  7  ,T  D ,  TP  T,7 1 0  0  ,YIT) - - - 

PRINT  4?,  <1,  7OIM.T0  ,IPT,T10 

-FauTUHi,-ro;or  pts~»  m, 2y.~0R0ER«-^Tt7; - 

C  2<*  “I  Ml  T  TIRE*  ".G12.5//1X,  ”RE  AO  Et/EPV  ~,I2,“  °  OINT  “» “»  Yl 

C  •  51  •  -  ■—  - 

K7*TO/0,03 

k*<i-<t  - — - - 

OtL*3»3T*IPT 


2(6,7) 
0)  — 


-wr 


0*  512 


HE1DI1, .3) (I  , COM4, A  Z, AZD, 0UM1, EL, ELD, EL  CO, 3 JM2, ELHN,  DUH3,  •  LSD, 
- C  1*1,  XT)  . . . .  .  -  . . 

uC  ♦  4  1*1, K 

- REACH.  ,3 1  T  ,CUH4,  A  Z,  AZO,  0UN1  »1TLiEX3,  ELU0V3'JNZyr  LMN,  DUM3, E  L3  D 

IF(EJc,U)  4  4,  46 

TTJ - FORMAT  (I2G12.S) - - 

45  lF(*n;  ( I-  1  ,  X  FT  >  .NE.  0)  GO  TO  44 

Ih*  (I-U7IPT  *1  - - - - - - 

IF(  H.N-t.  t>  GO  TO  49 

A  Z3  0  *  A  Z  ) 

TT9 - CONTINUE - 

A  T3<IH1*0JM1 

. X  TIM)  *3  JM2  •  •  - - —  - -  ' 

till Im(*£l-ELMN 

- —  —  S  (TH)  *  3  393  -  - - - - — - - - - -  •  • 

nA'XIP)  *7.5 

- XFr3UM4.LT;2877Vt-9AN(IM)*OOM47-(930-r-Tl-9^0Ji4) - - 

.  M*7  AN ( I (3 / jE  L 

, - -  XFtn.LE. Ml  GO  TO  44 - — - 

i  SET  *  I  SET  *1 

IF(I3£T,N£.l )  GC  TO  44  -  - 

AZT«  l*  A! 

.  VI - CONTINUE - - - 

R  H3  (  2,46  3  (T  IME.JUMS.TEL,  0UM7 ,T ELSO, 1 0*  l.KT) 

CC  A’  1*1, K  - . - . 

KEV:  (2,  .61  T  IME ,  GUM5,  TEl,  0UM7,TELSD 

-  -  -  IF(CTc  (?)  )  46,46  -  ” 

40  F  C*MA  T  (  >S  12  «  51 


63 


- RSnsiTiT-ri - 

htlO* ,  ULPHA  <I>,I*1,NPM1 

_  PHNT»,JPS,  <ALPM<II.J«1.J(PAR) _ 

*Aj *070  —  -  -  - 

ID  OC  8  1*1,  NPAR 

tVTT*.  l 

£(II«0.3 
A  (U  *  2 Vi 
COATINU; 
fT*SU8lT*tf 

V  29*4ZTAU-(A20-»100I 

ULTI Sit  G  aCPhST*  JT - 

C11J*AJ 

print  4 i'.talthjt ( np tts p* i •  npart 

FCRNAT  (“OALPHA*  G 12 .5,  “  l“,  G  12.5,  "I“,  S  (2 . 5,  “  I  ",  G12.  5.  “I”, 

"C  Cl 2. $ .'« “,G  1275,-1  G12.5)  ' 

PRINT  41,  AJ,IT,SU6IT 

TMHur-"  T«  ‘T(iTZ7S7"~n  our  iiiHSi_-7r5rwsTOin srrrr  ots«~-,t5i 

00  100  <«1,NPAR 
T03IT»1J9IITI 
OC  1?  L*1,NPAR 

i  k  to  «al  Fnn  cr  ♦ron*Jrnr,CT - - 

IM(L.G-:.2».AN0.a.LE.4>>  GO  70  12 

imLPRTN nrr.a .» roj  alpwa ntli  — althihiu - 

12  tCSTINU- 

C  RuTlN  TG  (AIFKANTAJ  I 
If  (AJ.JT.OLCJ)  GO  TO  20 

- Oil  A  j  - ' - - - - 

PRINT  43,  CAL  PhAN(NP )  ,  NP*1 ,  NPARI 

FRr*iT'4t;  aj,  it.sustt 

£  l<)  *0(0  ♦6«K> 

- E  «  >  *  S»  “  <  K  i - 

0C  15  M*l,NPAfc 

- SIP THTfR*  •  ALP  W UTTO 

15  CONTINUE 

- irrr«i  ;gt:t;5t  iKim - 

6C  TO  25 

73  ETO«^75»EOrT  ~~ 

IF  (MO.LE.1.51  A(«l*0.0 
75  Ck*B70‘ 

C.C  31  L*1,NPAR 

'  re  ltd  i  ’  r.  w.cien— nr~»rB - 


A 


..  40 


41 

11 


ir 


Cwi.O 

'TCRirsar 


IF  HK.  IE.t.  01  GO  TO  ISO 


4Uii*Tu«i1f7T 


00  3?  N*1,NFAR 


TcrrrsTnr?*  i  ti  ,m^*r 

SL12*$UN2*XI  (2,  N) **  2 

■jwrrhtfi - 


-jr 
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X1*S1RT(SUN1) 
"x2»51RT(SUHZr 
X3*Xl/X» 


TftIRT~3 1,  OLD  J.'UCPffATHPr,  PP«T,NT»AT7)",Tr;Tr- 


FCINAT  C»1J»  *,G12.5/»  ALPHA«(*,G12.5,*I»,J12.5,*I*,G12.5 

~I*  i*7Sl?Y5,  *  1*7617.  !.«i'fa2.VT*  ,  G 12.  5  ,*»»/ - 

3*  <Ht>»  »,G12.S/»  XI  (1>/<I(  21*  *,G12.5> _ 


100 


6C  TO  110 
CONTINUE 


OC  TO  U 
S113=0.0 


■CC'flE'C«'17KWTir 

OC  115  1*1*  NFAR 

TUTitmr- - 

F2(L,KJ*0.0 


TT5 — cxvrrrjr 

00  117  J*  1*  NFAR 

- DOT  1 7~  <•  1,'NFXR - 

117  X  It  J, 0  =  0.0 

- - cc"i»tr~r«i,~NFAR - 

120  SL13=SU13*(AeSt0(I))> 


~IF~(5ll  HT7CE  Sl'GO'TOTnrr 
OC  1  TO  1*  l.NPAR 
TC110_J«  1,1.  FAR 
Z(1«J)=0tN)*V(N«JI 

ciunoe - 

OC  1<*0  J*  1  *  NFAR 


lvO 


“DCT50  'Tilt  NFAR - 

OC  lbO  <* J, NFAR 
“XTTJ7  n  =  XI  U  ,1}  *?  TK  , XT 
continue 

LI  C  l  0.5 - 

OC  150  J*l,  NFAR 


.  ttd — sut  lv  nnc  nrrrrrji^*? 1 

,  SUHS31T  (SU7*) 

, - IiCT5F~J*T»  N  FA  R - 

.  155  *  (1,11=0(1,  J1/SUH4 


159 


lbO 


T7TT 


“RC'JNT  *2 
b*KO'JNT-l 

-utrrrrx-Tro - 

OC  150  0*1, NFAR 

'Fni.nii  ♦  - xr < koont ,o *7 txrcr 

C04TINJ5 

“DC  170  1*TVWFAR - 

F2(3,M)*F  (Q)*V(K,H)  ♦  F2ta.MI 

“ccrnmn - 

FH)  =  0.a 

'CCITIHUr - 

OC  190  1*1, NFAR 

“B  «  0  J  NT ,  r  )  *  XTOCCONTiTI  -F2  (Q,  IT 
SU15( 0) *0 ■ 0 

~0GT3(P’1*TTNFSR - 

SUIStQI  *SUM5  IQI*B(K0UNT,H)*»2 

- scmtDMsaRT  tsuHsnri - 

to  215  1*1,  NFAR 

- V  KOUNT.M)  *£.  (KQJNT ,  NT/S0N5Tm 

215  CC9TINU* 


150 


2M  0 


'KCiNTAOtnrrn - 

IF  <<3UtT.Lt  .NFAR)  GO  TO  159 

^*IT*1  - 

itniT*  j 

‘  FK*  1 - - - 

CO  250  <*  1,  NFAR 


TKT*n— 
OK1=C.O 
"a  rxT*2.r 


1 


I  •  t(  lb)  .  £&  ( lb)it61NT  (16).  EM  h,  41  ,EAI  NT  <4  ,4),  r  (  4) ,  CQC<  4,4),  C  I  4) 
2.XM10  0J)  |XM10  801  •  EOH(10  001  .EOOHUOOOI 
tG'lIXAL^NCE  iaflil)  ,  <EAU,1),EB  «nr.'<“AlMT  (1.1 J.E6INT  (1)1 

oat*  wr/i./ 


INITIALIZATION 


GE*COS  (ELG(KK)  J 
aa*.*-c£l»Alpha‘(  i)V?S 
IF(**3.GT.-200.I  S1»EXP(4RG1 
XJ(<l)*K3(KK)*cD0<<l<>*0£L 
1F(*L»H4(1>  .EG.  0.1  GO  TO  4 

i(4<  iriTis  i*  iJcnoiK  i  ieoofwam - 

(.C  to  1 

X 4(Kl)  =>XL  (  Kif  r*EG0 (XK)  *0£L 
continue 

t  c-n  <i)3  xn  or-x  nx  i  • 

IF(<<.u£.2)  £  CjH  (KKI  » C  EOH  (KK)“EDH(KK-1)1  /OEL 

"~K=T#T(K<)70EE - 

IF(«.LE.M)  GO  TO  10 

“rsjTsrsim — - - 

IF(I3£T.NE.l >  GO  TO  10 

Truiir  13'  CS  “  - 

h  (»)«T20*Cfl 


P<»,?)3(0.0071S«7*C8)»»2 

hs>*«  ^  mi  iv  c&-  xekp  cm » rz*nm  j/cB-nrnn »» z 

SSOMSORT  (P<1,  ln/CB-SEHPIll)**- 2*  (SORT!  PC  2,  21 1 /CB-TS  Cl>»**2 

T5r«x<  '  .  . . 

CCUINUi _ 


T6«-TnEI0»TAMSLGCKK)  I 
SC«LMC0»Ctfl*‘  2/OEL 

T«T  0;E - 

42»l.-<rAmKH-TAU<«X>)/3EL 
~C1R«C0»AI  PH*  CZJ  *CB 
A  (It  II  ■-COR»Tb 


A1»,N01I«-C0A*A2 

- *T!7NrjiA'nTZ)»»2 - 

A  It,  2)»T9 

x - 

C  COMPUTE  TB  A  IS  IT  I  CN  MATRIX  EA  AND  ITS  INTEGRAL  HINT 

*C - 

CALL  OSuRT  (NCIMtb«9ELtEB«  EBINT,  5) 

- CGW<r*ACT»H*  (4T»CB  - - 

C 

"C~ST  ART  "  I  NT  I G  RA  T 10  VCUOP - 

C 


C  1  ST»T(P(lt  1) I.SSO 

"99 - FtRHAUSGTZ."5) - “ 

C 


DO  120  J*  1»  NCIN 

- DTTlitMTT  *cA  (mmrm - 

120  CONTINUE 

1X0 - CONTINUE - - 

F  (1)  *  IC3>CR4)*X3(KK)  ♦CR4*X4( KKI 

- rt?T*x3(ificrns - 

IF(<*.GT.MI  F(  2)  bF(  21  ♦CR4»  <  X*.  (KK-M)  -X3<  XK-N) )  *  A? 

- 00*130-  I«irNTIH -  - 

UC  140  J»l,2 

- Dfll^om  *EA INTU,3r»FrjJ - 

140  CONTINUE 

- r  in«  cm - — - - 

C(IMO. 

TTO - C  CNU  NUE - - 

c 


C~  CONFUTE' TTCV 1PTA NOT  HATE IX - - - 

C 

- CCIITI)*  T*CFHriST*AtT>HAT«n,-AB3TnWtKKn»AL-PH«  tT  T*AJS  (EODH ( <K)  )  I 

1  •ECAl 

- CCrT2t2l«ALPHA(5l»SCAL»A2>»2 - 

IF(<<.GT.M)  COC  (2»  2  I*  I  ALPHA  (5)  ♦  ALPHA  <  6)  *ABS  C  EDHf  KK-M))  »ALPHA(7>  • 

- 1  '1ST  ED3H IKK-HI J  J*SCAL*A2**? - 

CALL  MULT  (EA INI  *CQ0  «N0IM(  N1,P1,  10) 

- CALr-PUCT(EAtPfKOIN,Nl,»2,10l - —  -  ' 

DC  220  I*  lv  NCI N 

- CC  -220  ■)«  I  •  KCIH - 

Pit,  I  )«PHI,  J)*P2«I*J) 

7ZV—  continue  -  - - - - — - 


>-TS  tCl-IST))- 


O-'COHPUTE  ERR3P  DUETff  ST ANOARO  DEVIATION 
C 


KTO  CONTINUE - - - 

E  J»  DEL*  (SHEA  N*KT*SSO) 

- RETURN  . . 

EM 
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iU5*0jTIN£  hULT(£,F,L,Ll,  M,MR) 

0  HE  MS  I  IN  E.  <  1)  » F  1 1)  iGiniiHIl) 

OCTO  r«i»L 

11*1  _ 

OClO  K*ltL 
UiF*0. 

nnro*i7CiTE - 

1clP*T£IP»t(J)*F<II> 

i  i»  r  t  *  i 

KK*(<-1)»L»I 
h(«»  *T;MP 
lb  G  €<< )  *T£MP 

ZnM%.n.llkETU4N 
LC  20  1*1 tL 
CC  20  <*I,L 
TE1P*0. 

II**' 

00  15  J*I  «L1  ft 

- t  rin»TFip  »&  i  jr*mn - 

lb _  II*TIfL  _ _ 

2b  h«K»*T£MP 

i.  2*1-1 . 

OC  31  I *1 »  L  2 
1 3*  I  *  1 

CO  50  J*L5,t 

ki« 

K2=  (  J*  1)  •  L»I 
30  H  « H  *  H  (K  21 

tn 

- bXfsoun  n  nrairr  i  N3ih  ,  atd  ewea.  eatnttnt  i - 

„  0HENSI3N  A  (l)tEA(lltEAINT(l)  tC  CEF  <  30) 

1C  S£is  EA«tXMA»0£O  fEAlNLlNTEGRAL  EA~ITT3  TJEC 

„  N  Cl  PI*  N  31  N*1 

.  NMN5H»N0IH 

.  NH1*NT-1 

.  CCiPfNT»*l. 

CC  13  I* 1  *  NT  PI 
«  II*  NT -I 

-  10  CCE7  ( 1 1 I *OEL*COEF<I  I»l)FFLOATII) 

«  "C  ' - _"“NT-«UST'eE"»T~ie<Sr“3 

„  CALI  0UG(NDIM,tAINT,A,:0EF(l>,COEF<2>> 

.  OC  68  t*3»NT 

CALL  KUlT  (A,tAXNT,NOIM,NN*EA,  1) 

«  lFCL.tQ.NTiGTT  T<T70 

..  bO  CAL  CI1G(NUXM,EAINT,EA,1.0,COEF<U) 

..  70  OcTMT  It*  1 V N Kf  NOIM1 

£A(II)*£A(II)*1.0 

50“  cnmxrj 

..  ENT 

SLT<rouYtNE~"oIt<n  NOTmi  *TB7TT7C"2J 
0 IiE  NS  I  3N  A  (  1)  t  B  ( 1) 

..  ltd*!*  Nil  H*1 

NN*N1IH»NOIM 

« - NRX^T-f-l - 

.  11*1 

.  iFcerreirrm  ~so  7o~ro 

.  _ CO  5  J*lf_NN.  N0IJ1 _ 

.  N*J*JMi 

,  CC  b  I  *  Jt  K 

- j—k  (ITirI  lgTn - 

.  A  III)  *ACII)  *C2 

.  ^  IX*  I  l*Nll  Hi 

. _ *ETIHN _ 

10  OC  r  JaltNNfNOXN 
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APPENDIX  C 

LISTING  OF  AN  AAA  GUNNER 
MODEL  SIMULATION  PROGRAM 


M8H,T20,CH70000.  L760!95,HEI, 2553960 
COMMENT. *NEM0WLSTMU6, 1 0*L 760295,  GY* 1* 

CO HMENT.  *AAA  NOOE6  3LANKINS  5IMJLATI0N  PROGRAM* 

AT  tACH.TAPEl,  0WL6MELSUB  J33, 10*L7  50  295,  CY*1 *NR* 1. 

FT  N, 

LGO. 

PROGRAM  S IM  J6<  IRPJT, OUTPJT, T APE1 » 

COMMON/S/CO  (2)  ,DEL  <  IH,NOtM,  T10  (21, X 3 (21,  EL,  EL  03,  AZOO,MTAU,RA, A2.NO 
Al,N02,NAA,USL.UAZ.ELTR,AZrR,ISET  , Zl,  Z2,  T  AU,T S  (1  5)  ,  T  E(  15)  ,T, IBL 
C  " 

C  THE  PURPOSE  OF  THIS  PROGRAM  IS  TO  SIMULATE  AN  ELEVTN  AND  AZ INUTH  TRApK 
C  TASK  IN  THE  TRA  IE  P-3 I°ECTE3  FIRE  (NODE  6)  SYSTEM 
C  SUBJECT  TO  OPTICAL  BLANKING 

C  INPUT*  THE  ELEVATION  ( EL)  l  AZIMUTH  I AZ)  ANGULAR  ACCELERATION  OF 

C  TARGET.  AND  BLANKING  OURATIONSIUP  T3  15)  IN  CHRONOGICAL _ 

C  ORDER 

C  OUTPUT*  ME4N  ANO  STAND  OEV  3?  LAG  ANGLE  _ 

C  *****  ALL  ANGLES  ARE  IN  UNITS  3F  RADIAN  ***** 

C  TAU*  DELAY  IN  SECONDS 
C  ALPHA*  PARAMETER  VECTOR 

C  ELERR*  MEAN  EL  LAG  ANGLE  (I.E.  TARGET  ANGLE-BARREL  ANGLE) 

C  AZERR*  MEAN  AZ  LAG  ANGLE 

C  ELSO*  STANOARO  DEVIATION  0*  ELEVATION  LAG  ANGLE _ 

C  AZSOt  STANOARO  DEVIATION  OF  AZ  LAG  ANGLE 

C  ELTR*  MEAN  EL  TRACER  ERROR  (TARGET  ANGLE-TRACER  ENOING  ANGLE) 

C  AZTRI  MEAN  AZ  TRACER  ERROR 

C  ELBAR*  MEAN  EL  DARREL  ANGLE _ _ _ 

C  OCL  *  TIME  STEP  USED  IN  THE  INTEGRATION  ROUTINE 
C  TS<  I)  *  STARTING  TINE  OF  I-TH  BLANKING  DURATION 
C  TE< I) I  ENOING  TIME  IF  I-TH  BLANKING  OURATION 
C  Y10I1 )*  INITIAL  GUESS  OF  EL  LAG  ANG.E  _ 

C  Y102*  INITIAL  GUESS  OF  A?  LAG  ANGLE 

C  UEL*  EL  CONTROL _ 

C  UAZ*  AZ  CONTROL 

C  COUI*  EL  RATE  CONTROL  COEFF  _  _ 

C  CO*  2)  I  AZ  RATE  CONTROL  COEFF 

C  Kl*  NO  OF  POINTS  IN  HE  ENTIRE  TRAJECTORY  _ 

C  K*  NO  OF  POINTS  AFTER  THE  first  TRACER  ROUNO  IS  FIRED 

C  ELOO*  EL  ANGULAR  ACCELERATION  3F  TARGET _ 

C  AZOOI  AZ  ANGULAR  ACCELERATION  DE  TARGET 
C  X3t 1) 1  EL  ANGULAR  VELOCITY  OF  TARGET 
C  X3I2M  AZ  ANGULAR  VELOCITY  OF  TARGET 

C  Xh*  ESTIMATION  ERPDR  OF  ANGULAR  VELOCITY  OF  TARGET _ _ 

C  EL*  EL  ANGULAR  “OSITION  OF  TARGET 

C  H(l)«  MODEL  PRE3ICTE0  LAG  ANGLE _ 

C  Ml  2) I  NOOEL  PREDICTED  TRACER  ERROR 
C  Pit, II  *  VARIANCE  OF  PPE3ICTS0  LAG  ANGLE 
C  P(2,2) I  VARIANCE  0*  PREOICTE3  TRACER  ERROR 
C  TO*  THE  INITIAL  FIRING  TIME 
C 

REA3».K1,T0.IPT,  IBL _ 

PRINT  3,K1,T0,I3T 

3  FORMAT ( LM1, -NO  3F  PTS*  ", t *, 2X,*INIT  TIME*  " ,G12.5//1X,"REA0  EVERY 
C",I2.“  POINT"/) 

IFdBL.GT.O)  READ*,  ITSIO,  TE(K),K*1,  I9LI 
PRINT  ll.IRL 

11 _ FORMAT  (IX,  I5j  IX j  "SLANKING  INTERVALS  ABE  "/) 

IF(  I9L.GE.  1)  PRINT  i»,  JTSIO  ,  TE IX)  ,K»1,  IBU 
«,  FORMAT  I5UX,"(",  F9. 2,",  ",  F  9.  Z,  "I  "I  I 
KT-T0/0.03 
K»K1-KT 
T*T0 
N01M**» 
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X  PRINT  *20 /t®T 
C3(l)*1.3<* 

CO (2) =1, 2  6 
ND1*N0IM-1 
N02*N0IM-2 

NAA=N0IH/2-l _ 

rsr=i 

DEL=0.03»I°T 
ELSO*.005  *•1.5 
AZSO*.OQ5*»1.5 

zi=o. 

_Z2*0_, _ 

IS£T*0 

UEL*0. 

UAZ=0. 

PRINT  7 

FORMAT  (/1H  ,2X,"TIME'*»1X»'*EL  YE.  “  ,  9X  .  "EL  ERR  5  X,  “E  LSD" ,  6X, 

1**EL  CTR**.  6X,  **A  Z  VEL‘‘«6X.‘*AZERR**«6X,~AZS0‘*.  6X,’*AZ  CTR" _ 

2,6X,**EL  TR",6X,“A7  TR*V) 

REAO (1,2)  (T.OUHt ,AZ,AZ3,AZOO,£L,£LO,ELOO.AZMN, X  ,ZSO  ,  S,  I*  1,  KT 
C  > 

00  5  1=1,  K 

REAOU.ZIT.lUMl.AZ.AZO.AZlI.El.iLD.ilDD.AZMN.X.ZSD.S 

IF (EOF (11)1,1 _ 

IF(H00(I-1,IPT).NE.D>  33  T 3  5 
IM® (1-11/ I°T*1 
T*T0*( IN-1) *OtL 
TAU=7.5 

IF(0UM1.L£.2877. )  TAU=(1UN1/ (  930.  - .  13*0UH  II 

NT  AU®TAU/0EL _ 

IF (1ST . EQ< II  DTA  U*  TAU 
IST®IST  *1  _  _ 

RA*NAA/TA0 

A2=l»-(  TAU-9TAUI  /OEL 

OT  AU*TAU 

X3(  1)*EL0 _ 

X3 ( 2 I  *  AZO 

IF((IM-1I  .Si.NTAUI  G3  TJ  3 
TAUR»TAU 

IFIQUM1 .LE.  ALOII  TAUR*TAU*ANAXl(0.6,OUNl/5dJ0.l  _ 

Y10I1I  *-TAUR*EL3-.  001*(5.2»TAUR*.<.85*TAUR**2I*COS(  EL*0.05I 

V10  2»-0.0  25»SI3S(1..AZQ) _ 

Y102*-TAUR*A73 

IF(  IH, N£,  1)  GO  TO  13_  _ _ _ _ 

ELO®EL 

AZO»AZ  _ _  _ _ 

E10*Y10<1) 

_ 

GO  TO  10 

ISET-ISETM  _ 

IFdSET.NE.il  33  TO  10 

Z1®EL~  (EL  O-ilOl  ♦.901*(5.2*rAU*.L86*TAU**2)  *SOS(  ELO  ♦  0^051  _ _ 

Z2*AZ-(AZ0-E20) 

_  _C  ALk_QJ SSEiALiaSRiitSOI _ 

ELBAR«EL-ELERR 

IF(IIM-l)  .Li.HTAU)  Yt0(2)«flD2*3OS(IL6AR) 

CALL  03SAZMA7ERR,  AZS0.EL3ARI 

IF((MOO(IH-t,I»?IYT».E7.0>  .3R.IIH.E3.il)  PRINT  6,T  ,  X3( y_t_ELERR. 
l,ELSD,UELtX1(2),AZERR,AZS0,UAZ,ELTR,AZTR 

FORBAT  ( 11G1Z.5I _ 

CONTINUE 

FORMAT  (12G12. 5)  _  _  _ 

STOP 

CNO 


SUBROUTINE  13SE-6IELERR.ELS0I 

C0MN0N/SIC9  (2)  ,1EL.<K.90,V10f 2I.X30(2I  ,  EL.  EDO,  A  ZOO  ,H,R»,  A  2^ 


o  o  o 


A,N01,N02,NAA,U,UA?  ,  ELTP,  »ZrR,ISrT,Zl,22,TAU,  T  S  ( 15  ) ,  TE 1 15  1  ,  T  ,  IBL 
OINENSION  W(4),P(4,4)«Pli4,6)»P2(4,4),ALPHA(  7)  ,  A(  4,  4) 

1, B(16)  ,EB(15)«E3XVTtl6),£A(4,4),EAI'tT(4,4),F(4)  ,CTC  (4,4)  ,  C(4) 

2,  X3(  1000)  ,X4I132  0)  ,EDH( 1003) , ED3H( 10 00) , THET  (130  0)  ,ALP5 0(1000) 
EQUIVALENCE  ( A  ( 1 , 1  >  ,3(1) )  •  (E  A(l,  1)  ,  •  B(  1)  ),  (E  AINT11 , 1)  ,EB€M  ( 1)  ) 

DATA  NT/It/ _ 

DATA  ALPHA/1. 547 1,  .  017491,  .  024433,  .42318,  .22 446E-7,  .  1797SE-J, 

A  .17302E-3/ _ ; _ 

INITIALIZATION  .  _ _ _ 

if((kk-d.g;.h)  go  to  3 _ 

if(kk.gt.ii  go  ro  5 

N1*N0**2  _ 

NOIH-NO 

ALP1*ALP10*ALP-Ii  (l»_ _ _ __  _  __ 

ALP2*ALPHAI2) 

ALP3*ALPHA<3) _ 

ALP4«ALPH  A(  4) 

ALP5*ALPHA(5)  _  _ 

ALP6«ALPHA(S) 

ALP7*ALPHA(7>  ___  __ _ _  _ _ _ 

ALPS6(1>*AL*5 

ACC-0. _ 

IFLAG*0 

SCAL»Ctl(ll  **2/3EL _ _ _ _ 

00  1  1*1, NO 

00  1  J*l,  NO  _ _ _ _ 

P(I, J>*0. 

CQC(  I,  J)»6. _ 

I  A  1 1, J)*0 » 

00  11  1*1, NO  _  _ _  _ 

Htll*0. 

0(I)*0.  _ _  _ 

F( I) *0. 

II  P(I,  I)  *0.0030 _ 

P(1,1)*0. 0010296279 

P  (2, 21  *0.0300335677 _ _ _  _ _ _ 

X3(  1 )  *X30  ( 1) 

X4(  1)  *0.  _ _ _ 

EOH(l) *X3 (1) 

_ EOOH ( 1) *0  » _ 

S1«0. 

5  M(1)>V10(1)  _ _ _  _ _ 

6  IFdSET  .EQ.l)  41  2)  *21 

IF(IBL.LT.l)  53  TO  10  _  _ 

ALPS*  ALPHA  (5) 

IS«IFLAG»1 _ 

00  12  I*IS, I3L 

IF(T .GE.TS (I)«A'IO.T,LT.fs)I)l  G3  T3  IS 
.tt«ahini(i.s,a;c/3.) 

IF(T  .GE.TE  1 1)  •  A4D.T  »LT  .(TE(I)*ATTI)  GO  TO  16 
IFIT. GE.ITEII)*4TT>>  GO  TO  21 
_ ACC«0. 

IF(T.LT.TS(IS)»  GO  TO  18 
GO  TO  12 

16  ALPl«ALPie*UL<MA)l>-AL»l)(  •  ( 1  .-£<»( -C  .4  3*  (T-TE  ( I)  )  )> 

ALPS»-0. 0001* (l. -EXP (-0.43* (T-TE(IM I) 

GO  TO  16 

1 5 _ ACC*  ACC  *0  EL _ 

IFLAG*I-1 

ALP1*ALPHA(1I*E<P(-0.0*5»»(T-TS(I) )) 
ALP3»ALPHA(T)*E<P(-3.52*(T-TS(I> )) 
ALPS«O.OOOl*(l.-EXP(-9.12*(T-rS(II )) I 
ALP 10*ALP1 
<9  !9  46 . 


21 


AlPl«AlPHA(l) 

ALP3«AlPNACT) 

AC.P5*ALPHA<5) 

GO  TO  18 
12  CONTINUE 

18  CONTINUE _ _ 

ALP50C  KK)  *ALP5  ’ 

ANG«-0EL*AU»1 

IFCARG.GT.-200.I  S1*EXPCARJ> 

TMETC*«)=EL 
K1*KK*1 

K2«KK-1 _ 

C 

C  COMPUTE  TARGET  VELOIITV  ANO  ESTIMATION  ERROR 
C 

X3IK1)  sX3C<<)  ♦  EOO’OEC. 

ZFIALP1.EQ.9.)  30  TO  4 

_  X4CK1)  *Sl*X4<Kl<>  »;00»Cl.-51)/4i.»l _ 

60  TO  3 

4  X4(<1»«X4('<>0  *E00»0SL 

3  CONTINUE 

E0H(KKI*X3(<K)>X4I«)  _  _ 

IFCXK.GE.  21  EOON  C  KK ) *  C EON(  <K> -EOHC K2 > ) / DEL 

_ X30  Cl)  «X3  CXI) _ 

C0R*C0C1)*ALP2 
C  R4*C0 ( 1) • AL°4 
IFI CKK-1)  .L5.M)  GO  TO  150 

All  »1 .  *0 •  0  0 1*  I  5,  2*TAU*0.48jylTAU»*2)  »SINC_THET  UX^M)  ) 

A(1,1)*<C0R 

AC1,2)»-G0(1>»4LP3 _ 

A  ( 2, NOl) 3 A  C t,l)*ALl*A2 
AC2»ND)*A(1,2)*AL1*A2 
00  2  1*1, NO? 

J1*I  *2  _  __  _  _ _ 

ACJ1,I)*RA 

ACJ1,J1)*~R4 _ 

2  CONTINUE 

CALL  OSCRTCNO,9,OEL,EB,EBI'IT,5)  _  _ _ 

CR3*C«4*AL1*»2 

SCAL13SCAL*CALI*A2)**2  _ 

C 

C  COMPUTE  MEAN  TRACKING  ERROR  CI.E.  LAG  ANGLE) _ 

C 

U*AIP2*HC1)  ♦ALB3*tCCZ)*ALP4*  C  X3  4JJCK ) ) 

00  110  1*1. NO 

00  120  J*1,N0  _  _ _ _ _ 

OCI)*OCI)  ♦EAC  I,  J)  *NC  J) 

120  _ CONTINUE _ 

110  CONTINUE 

F<ll*Cl.-CR4)*X3<i<KM-CR4*X4CKK> 

FC2)*X3CKK)  *CR3*  CXACXK-N)  - <3  C  KK-M)  )  * 0 . 00  1*  ( 1.- A 2)  •  C  5. 2*  0.  972*f  AU 

1)*C0SCTHETCXK-N)  )  _  _ 

00  130  1*1, NO 

_ DO  140  J*  1,2 _ 

OC li *0(1) ♦EAINT  CI,J)*FCJ) 

140  CONTINUE 
MCI)*OCI> 

OCI)»0.  _  _ _ 

130  CONTINUE 

c _ 

C  COMPUTE  COVARIANCE  MATRIX 

C  _ 

CQC  Cl,  D*  CAt°5*At*6*ABS<  EJVC  KK) )  ♦AU*7*ABSCE00HC  KX)  )  ) 

1  *SCAL 

COC  (2,2  )*  CALPFOC  K<-M>  ♦AL»5»A3S  CEO-C  CXK-M)  )  *»L  P7»  ABS  C  EOOHC  KK-H) 

1»  *»SCALi 


call  MULrtEiiNr.coc.NO.'u,  »i,ioi 
CALL  MUlTCEI.P.MO.NI.PI.IJ) 

00  22u  1*1, ND 
00  223  J*t,V0 
PII,  J)*P1CI,J)*32(I,  J» 

220_  CONTINUE _ 

150  CONTINUE 

ELERR=HCi)  __ 

ELSD*SQRT C»C1,1) > 

ELTR*HC2) 

RETURN 

_ ENO _ 

SU3R0UTINE  38SAT6C  AZERR.AZ50, ELJ) 

COHNOM/S/’CO(2I  ,3EL,«,ND,U0C2>,X30C2),EL,ELOO,A2OO,H,RA, 

1  A2,N01,N02,NA1,UEL,U,ELTR,AZTR,  ISET  ,  21,  Z2,  TAU,  TS(  1  51  ,TE  ( 15)  ,T,IBL 

DIMENSION  NIL)  ,’(<»,<*)  »°1 C  4, 4)  ,P2  14,4)  .ALPHA)?),  A(4, 4)  _ 

1,B(  16)  ,£EC15),£3I'|T(16),E4C4,4>,£AIVTC4«4),FC4>  ,COC  (A, 4)  ,  0(4) 

_  2«X3<  1000)  .X4C1H0I  .EQHC1003)  .ED3H1  10  00)  ,  ALP50) 1  BOO ) _ 

EQUIVALENCE  C »  C 1 , 1 )  ,  3)  1) )  ,  (  EA  C 1, 1)  ,  EB  C 1)  ) ,  (E  AINU  1 ,  1),EBIMC1)> 
OATA  MT/1 ,/ 

OAT  A  ALPHA  AS. 5 304,  .13*94,.  127  23, 1 . 03 5 3, . 2528 6E -5,  .  19766E-3, . 757 85E 

A-3/  _  _  _ _ 

IFC(KK-l)  .GE.M)  GO  TO  6 

_  IMKK.6T.1)  SO  TO  S _ 

ni*no*»2 

C 

C  INITIALIZATION 
C 

00  1  1*1, NO 

_  00  1  J*1,N0 _ 

p<i,  j)*‘o. 

CQC  (I,  J)  =  0  . 

I  A (I, J) *0, 

00  11  1*1 ,N1 
M 1 1) *0, 

O) I)*0, 

>ri)»o, 

II  PCI,  11*0.0190 

P) 1,1>*C0.'JJ59183*C3SCEL3)  I  *»2 
PC2,2)*C0.31T155r*COS(£Li))**2 
X3  C 1 )  *X30  C2) 

X4C1)*0. _ 

EON ( 1)  *X3 C 1) 

EOOHCl)=a. 

S1*0. 

ALPt*ALP10*ALf>HA(l) 

ALP2*  ALPHA)  2) 

ALP3*ALPHA(  C) _ 

ALP4*ALPHA(,I 

ALP5*ALPHA<5) 

ALP6*ALPHA(5) 

ALP7*ALPHAC7> 

ALP50C 1)  *  AL°5 

_ ACC«0. _ _ 

iPLAG*0 

C 

C  COMPUTE  ANO  STORE  STATES  X3  AND  XL 

C  _ _ 

5  Nll)*riOI2) 

6  _ IFCISET.EQ.l)  4)  2)  «Z2»S3S  C  EL6) _ 

CONTINUE 

C3*C0SCELG> 

IFCIBL.LT. 1)  S3  TO  18 
ALP5«ALPHAC3) 

IS*IFLAG*  1 
00.12  X*IS, IBL 
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^  wo  O  O  ft-  O  OO  ooo 


IF(T.GE.TS(I>  .HO.T  .UT.rEdll  33  T3  15 
ATT*ANIN1(1. 5,430/3.) 

IF<T.Ge.TEUt.4MD.T.LT.(TE(II*ATm  GO  TO  16 

IF<T.G£. <T£<T>  *ATT) >  GO  TO  21 

ACC*0. 

IF<  T »LT .TS( IS) )  GO  TO  18 _ 

GO  TO  12 

16  AlPl*ALP10»<AltMAU>-ALn3>Ml.-EXP<-0.a3MT-TE<I>  »  )) 

ALP5*-0.0  01*<l.-e*P(-0.43MT-TE< I>)>  > 

GO  TO  16 

15  ACC=ACC*DEL 

_ IFLAG»I-1 _ _ 

ALP1* ALPHA ( l)*E><P(-l.07;&*(T-TS<  I)  >  I 
ALP2=ALPHA<2)  *E<P|-3.12*<T«TS(II  I) 

ALPJ*ALP4A(T)  •  *<P(  -0 . 52*  f  T-TSt  I)  )  > 
ALP5*0.001*(1.-EXPJ-3.12*C  MStltlll 
ALP10=ALP1 

GO  TO  18 _ 

21  ALPliALPHA(l) 

ALP2* ALPH  A  ( 21 
ALP3=ALPHA<  Jl 
ALP5*ALPHA(5> 

GO  TO  18 

12 _ CONTINUE _ 

18  CONTINUE 

ALP50(KK>=ALP5 
A  RG=-OEL*AL3l *33 

IFIAVG.GT.-200.)  Sl*EXO(4Vi>  _  _ 

Kl- KK*1 

_ K2*KK-1 _ 

X3(K1)«X3(<>0  ♦4TDD*OEL 

IF (  ALP1.E0.3.I  OO  TO  A  _ 

X4<Kl)*$l*IX4<<10  *A200*3jLI 
GO  TO  3 

X4(ia)=X4(<'0  *4* PP*3EL 

_ CONTJNUE _ 

COMPUTE  ANO  STOVE  ESTIMATED  T4V3ET  _0EL03ITY  AND  ACCELERATION 
E0HJK1)=X3C<1) -TMK1) 

IF (<2.G£.  1)  E03HC<KI*(E3H(«)-E3H<K2))/OEL 

_ X30 (  2)  *X3<<1) _ 

IFUKK-l!  .LE.MI  30  TO  1S3 

OO  2  1*1 ,  N02  _  __  __ 

J1*I *2 

A  (  J  It  I)  *  VA  _ 

A  ( Jl  t  Jl)  *-V4 

CONTINUE  _ 

TmE80* (ELG-3EL 3) /DEL 
C9»COSIElG> 

Tb* -TM£BO*TAN ( ElG ) 

SCAL*(C0(  21  *3?l  **2/0EL 
CO  V*C0  <  2)  *  ,L°2*3B 

_  au.1)»-co»*ts _ 

A<l,2l*-:o«2) *ALP3*C3 
A 1 2t NOl )*-C)V*4? 

AI2,NO>*A(l, 21*42 
A)2,2)«TB 

COMPUTE  TRANSIT  TON  MATRIX  E>  MO  ITS  INTEGRAL  E4  INT _ 

CALL  OSCVTHO,3,OEL,E6,E3I'ir,5) 

CV4*C0«2»*ALPA»:8 

COMPUTE  MEAN  TRACKING  EVVOV 
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U*ALP2*W( 1) *AL»d.l ( 21+AIP4* IX3(KK) -X4IKK)) 

00  110  1*1, NO 
00  120  J^l.NO 
0<I>*0<I>  *EA<I,J)*W<  J) 

120  CONTINUE 

11 0 _ CONTINUE _ 

F«1>*<CB-CR4>*X3CKK)  *CR4*X4<KK> 

F(  21  *X3<  KK1  •CB*:R4*{XM«-1)  -X3<  KK-HI  *A2 
00  130  1*1, NO 
00  140  J*  1.2 

0(I>*0<I)  ♦E4INHI,  J)*F(JI 

140 _ CONT INUE _ 

wn>*b<i> 

0(I)*0. 

130  CONTINUE 

c 

C  COMPUTE  COVARIANOE  MATRIX 
2 

cocfi  ,T»*<al>5*al>6*abscs:h<kki  )  *a  l  >  z *  a  3  eTeoo  hTkkTTT 

1  *SCAL 

CQC  •2,2>*(ALP5H«-i>*4LP4»ABS«:0Mt<K-M)  ) 

1  ♦ALPT»ABS(EPOH ( < <-D ) I • S GAL » A2 •• 2 
CALL  MULT ( SAINT, COC, NO, Nt, >1,101 

_ C A LL  MULT (EA . P, N 0, Nl, P2, 13  1 _ 

00  220  1*1, NO 
00  220  J*1,N0 
P(I,J)*P1(I,J)>>2II,J) 

220  CONTINUE 
150  CONTINUE 

_ OELG=ELG _ 

AZERRsMCU/OO  * 

AZSD*SQRT (°(1,1))/C9 

AZTR*W(2I/C3 

RETURN 

ENO 


DIMENSION  SI1),  C<1), 5(161, 411) 

00  10  1*1, L 

11*1 

00  10  K*1«L  _  _  _ 

TEMP*0. 

00  S  J»I,Lt,L _ 

TEHP*T  EMP  ♦SCJ)*e(III 
5  11*11*1 

KK*  <K-1)*L»I 
HCKK»*TENP 
10  G<KK>*TEMP 

_ IF<MR,EQ.1I  RETURN _ 

00  20  1*1 ,L 

00  20  K*I,L  _ 

TEMP*0. 

II*K 

00  15  J*I,L1,L 

_ TEMP;  TSMP*GlJ>»StII> _ 

15  II*II*L 

K«»<K-1I»L*I  _  _ 

20  HI <KI *TEMP 

L2*l-1 

00  30  1*1 «L2 

_ L3*J  ♦  l _ 

00  30  J*L3,L 
Kl*<I-ll*l*J 
K2*« J-1»*L*I 
30  H(K1)*H(K2) 

ENO 

SUfROlfTINE  OSCRT  (NOIM,A,OEL,EA,EAINT,NTI 
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DIMENSION  AC1),S  All)  ,£4IUI1I  ,C3EF< JO) 

C  SETS  EA*£Xs<A*OEO»SAI'(r»INT£GRAL  EA  0  TO  DEL 

NOIHl  =  NOIMn 
NNsNOIM*NOI>1 
NTM1*NT-1 

COEF<NTI*l. _ 

00  10  I=1,NTM1 
II*NT-I 

10  COEF C III *OEL»C3SF< 11*11 /*LD ATI  II 
C  NT  MUST  31  AT  LEAST  3 

CALL  OIAG(N3H,SAINT,A,:3EJ(ll ,;0£F(2)I 

DO  60  L-3.NT _ 

CALL  MULT  (AfEAIST  t  NO IM*  NN ■  £  A  *  1 ) 

IFIL.EO.NTISO  n  70 

61  CALL  OlAG(NOIM,£AINT,EA,t.3,C0EF<L)l 
70  00  80  II*l,NN,N3IMi 
E  Al  1 1)  =E  A  ( 1 1)  +1 .  0 

80  CONTINUE _ 

ENO 

SU3R0UTINE  OTASC NT  IN , A , 3 , 3 1, C2I 
DIMENSION  A  1 1 )  ,  3  (  1 ) 

N0IM1=N0IM*1 

NN*NOIM*NDIM 

_ N  MlaNOIM-1 _ 

II  *  1 

IF (C 1 • £Q< 1.3)  S3  TO  10 
00  5  J=1,NN,N0I'* 

K= J+NM1 
00  4  I=»J,K 

_ 4  A  I II *C1*3 ( II _ 

A(III*A»III*C2 

5  II=II*NQIM1 
RETURN 

10  00  7  J=1,NN,N0I1 
K*J*NMl 

00  6  I*J,K _ 

6  A(II*S(II 

A  I  II)  *A  I  III  *C2 

7  II=II*NDIM1 
RETURN 

ENO 
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